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Abstract - We show that Lyapunov exponent for the Sinai billiard is A = —2 log(i?) + 
> ■ C + 0{Rlog^ R) with C = 1 - 41og2 + 27/{2tt^) ■ C(3) where R is the radius of the 

I circular scatterer. We consider the disk-to-disk-map of the standard configuration 

■ where the disks is centered inside a unit square. 

o 

: 1 Introduction 

I ■ In this paper we will confirm, and make more precise a widely believed conjecture concerning 

O I the Lyapunov exponent of the Sinai billiard in the small scatterer limit. We will find that the 

2 ' Lyapunov exponent for the disk to disk map is A = — 2 log(_R) + C + 0{R\og^ R) and derive an 

I exact value of C for a particular configuration of the billiard, see below. The Lyapunov number 

is known to be well defined for this system . There has been a lot of good argument for the 
term — 21ogi? H, ^, but the constant term has shown harder to approach. 

■ The (largest) Lyapunov exponent is defined by 

1 

A= lim —log\A{xo,n)\ , (1) 

n — >oo Ji 

defined for almost all initial points xq. A{xo,n) is the largest eigenvalue of the Jacobian of the 
n'th iterate of the map. 

If the map is one-dimensional x i— > f{x) then 

1 1 "^^ 

A= lim -log|r(xo)|= lim - V log , (2) 

n — *oo 77, n — >oo Ji ^ — ^ 

i=0 

which we, assuming ergodicity, can rewrite as 

A = / \og\r{xMx)dx , (3) 
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where p{x) is the invariant density. The problem is that p(x) is usually unknown. 

For an area-preserving map, as we will study, the problem is the reverse. The invariant 
measure is known (it is uniform) but A{xo,n) does not have the simple multiplicative property. 

Our strategy will be the following. We will first define a multiplicative weight giving the 
correct R —> behaviour of the Lyapunov exponent. Then we will use the uniform measure to 
compute this limit. Central to the development is the distribution of recurrence times and the 
majority of the paper will be spent on its derivation. 

We will consider the standard configuration of the Sinai billiard (although some of our con- 
siderations apply to slightly more general cases), a unit square with a circular disk with radius 
R centered on its midpoint. We will study the disk to disk map, the Lyapunov exponent for the 
flow is available via Abramov's formula Istl. 



2 Calculations 

2.1 Reducing the problem to one iterate of the map 

We are interested in computing the largest eigenvalue of the Jacobian of an arbitrary long tra- 
jectory. For the time being, we consider a rather general dispersive billiard with a circular disk 
with radius R inside a polygon and the associated disk to disk map. The disk does not touch the 
polygon. The Jacobian associated with one iterate of the map is, see e.g. 




1 -I 'Hi 

i?. cosav * 
2 ^ 

R cos cti 



(4) 



using a coordinate system following the trajectory - y being the transverse coordinate, li is the 
traveling distance and is the angle with respect to the normal. We have omitted an overall sign 
of J which depend on the number of times the particle hits the polygon during the disk to disk 
flight. For an area preserving mapping it is sufficient to know the trace of the Jacobian in order 
to compute the eigenvalues so we are interested in computing the trace of a product of matrices 
n"=i Ji- We want to expand this trace in power series in R since we want to extract the R ^ 
limit of the Lyapunov exponent. To this end it is convenient to split up Ji into the following sum 

_ 2h ( f 1 \ ficosa, / 1/k I \\ 
^ RcosaAy l/k 0j+ 2 I V'. / J ' 



or 



with 



We write the trace as 



with 



2L f . Rcosa.i , . . 

J^ = -l^^\A + TT^B,} , 6 

Rcosai [ 2 ' 



i=l i=l ^ ^' 

n „ 

^n = Trl[{A + ^^B,) . (9) 
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The idea is now to show that is bound close to unity when R ^ and n co. First we use 
the fact that all matrix elements are strictly positive and cos a < 1, so 

71 n „ 

l = rr[]Ai<*„<rr[](A + -Bi)=*n • (10) 

i=l i=l 

We expand in powers of R 

i?" , (11) 

where for instance 

ao,n = 1 (12) 

and 

ai,„ = Tr(Si^2 . . . ^n) + Tr(AiB2 ... ^„) ... + Tr(Ai . . . (13) 

etc. The coefficient « is a sum over ( ) combinations of A's and B's. A trace of such 

m 

a combination in a^.n is a sum of terms like . . .li^). As all U > Imin such a term is 

l/ik^k^ ...lij}< 1/Cm- We can thus write 

R R R 

$n<l+ai,n(^ ) + am,n{-, )™ + <n(^ )" , (14) 

''jnin ''min ''min 



where 
with 



oi,„ = Tr(BA"-i)+Tr(ABA"-2)... + rr(y4"-iB) , (15) 



Let us now take the trace of a particular combination of A's and 5's. This can always be organized 
as 

rr(i^iB'=M-''^B*^= ...i-'^E'^^) , (17) 

Y.a=li.ja+ka) = n , Ea=l = ■ (18) 

The matrix A obey = A so we reorganize the trace above as 

Tr{AB^^AB''^ ...AB'^'-A) . (19) 
We also have AB'^A = [k + 1)A so the trace above is bounded by 

L 

Tr{AB>'' Ab'^K . . Ab''^ A) = JJika + I) < 2"^ . (20) 

a=l 

So we finally have the following bound on <!>„ 

1 1 +( ; )(y^) + ( " ){^r +•••(" ){^r = (1 + 7^)" • (21) 

''min ''min ''min ''min 

Let us now recall the definition of the Lyapunov exponent 

A= lim -log|A(a;o,n)| . (22) 

n— s-oo 77, 
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The trace TrJ(a;o,?^) for an area preserving mapping is TrJ(a;o,n) = A{xQ,n) + 1/ A(xo,n) so 
for large n we have, provided that the hmit above exist, |A(xo,n)| = |TrJ(xo,n)|+ exponentially 
small terms. Using eq. (0) we get 



12/ 1 
A= lim -logTT(-— ^— )+ lim - log 4- , (23) 

n^oo n RcOSai n^oc n 

i—l 

and according to eq. ( ^ we get < lim„^oc ^ log's* < 2R/lmi7i- and 

1 " 91 

A= lim -^log(-— ^) + 0(i?) , (24) 

n~*oo n i? COStt; 

z— 1 

saying that, to leading order, the Lyapunov exponent is just the arithmetic average of a scalar. 
This enables us to rewrite it as phase space average 

where ls(xs) is the travel length to the next disk bounce and oi(xs) the scattering angle as 
functions of the phase space point Xg on the disk. We choose as coordinates for Xs the scattering 
angle a and an angle measured along the rim of the disk. The invariant measure is then 
dxs = d<p d{sm a) / An . 

2.2 Relating the Lyapunov to the distribution of recurrence times 

We can now divide the phase space integrals above into several terms 

A = y \og{2/R)dxs ~ J \ogcosa{xs)dxs + J \ogls{xs)dxs + 0{R) . (26) 
The first term is trivially = log(2/i?). The second term is also easily computed 

rTr/2 



f 1 r'^ 

I log COS. a{xs)dxs — I log(cosQ;) c?(sinQ;) — log(2) — 1 
To evaluate the third term we insert the identity 1 = J S{1 — ls{xs))dl 



(27) 



^ loglsixs)dxs = J I dxslogils{x,))S{l - ls{x,)) = I dllog{l)p(l) , (28) 

where 

p{l) = / dxj{l~h{x,)) (29) 



is the distribution of recurrence times. We are going compute the small R limit of this function 
for the standard version of the Sinai billiard, with the disk placed at the center of a unit square. 
We will discover that if we formulate p{l) in terms of the rescaled variable ^ — 2RI then p{l{£,) 
will have a well defined limit (if we provide p{l) with some smearing) when R gets small, and the 
following function 

c(i?)= / \ogC-pm))d^/2R (30) 



will have a well defined limit and the Lyapunov exponent will be 

A = -21ogi?-log2 + l + c(i?) + 0(i?) . (31) 
The whole problem thus lies in computing the function p{l{S,)) in the limit of small R. 
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2.3 Calculation of the distribution of recurrence times 

It is much more easy to grasp the geometry of the bilhard if it is represented as a regular Lorentz 
gas on a square lattice. The midpoints of the disks are sitting on lattice points in an integer 
lattice. To each lattice point is associated a lattice vector q = (m, n) . An immediate observation 
is that only disks represented by coprime lattice vectors {gcd{'m, n) — 1 where gcd means greatest 
common divisor) may be reached. We send out a ray from the (0, 0) disk. The phase space of 
the surface of section can now be partitioned into subsets Qq] a phase space point Xs G riq if the 
trajectory starting at Xs hits disk q. 

For finite disk radii not all coprime disks are accessible, some ilq —&■ To deal with this 
problem we will first construct a scheme for enumerating the coprimes. 

2.3.1 Enumerating coprimes, or climbing the Farey tree 

Without any loss a generality we restrict the considerations to the lattice points (m, n) in the 
first octant, that is m > ti > 0. Because of the restriction gcd{m, n) — I there is a one-to-one 
correspondence between these coprime lattice vectors q = (m,n) and rational numbers n/m. 

The set of all coprime vectors q = (m, n) with m < M is called the Farey sequence of order M. 
Next define the cross-product between qi — {mi, ni) and q2 = {m2, ^2) as qi xq2 = min2~nim2. 
Consider a vector q — (m, n) an draw a line from (0, 0) to it. Call the lattice point closest to the 
line and above q', and below is q" the neighbors of q in the Farey sequence of order m. The Farey 
theorem then states that q x q' = +1 and (for symmetry reasons) q x q" = —1. Geometrically 
this means that the distance from q' and q" to the line is l/|q|. A consequence of Farcy's theorem 
and the fact that q is coprime is that q = q' + q". We will call q' and q" the mother and father 
of q. 

The coprime lattice vectors, or the rationals, may be organized in the Farey tree 0, see fig 
la. There is a link between two vectors q' and q" if and only if they are Farey neighbors, that is 
|q' X q"| = 1. From each vector there are two links going upwards (towards the root as the tree 
is drawn upside down), one up to the right to the mother and one up to the left to the father. 
One of the parents lie on an adjacent level, the other is not. 

We now claim the following: 

Every coprime lattice vector q, except (1,0), (2,1) and (1,1), can uniquely he written as 
q = q' -|- nqc where n > 2 and q' is the mother or father of qc 

Let's assume it is possible and show how to do the construction. Call q = q„ = q' + nqc- By 
assumption we have |q' x qd = 1. It follows that q„ is linked to qc for all n. It also follows that 
q„ is linked to q„_i which is linked to q„_2 etc. The genealogy must thus be as in fig lb. 

So, given a vector q we do as follows to uniquely determine n, qc and q'. First we identify 
the parents of q. One is on the level just above but the other is not. The rule n ^ 1 force us to 
choose the other as qc. Now we construct qn-i, qn-2 iteratively until we reach q„_m = qo = q' 
for some m. It is easy to identify we this happen, since qi is on a level below qc and qo is located 
above it. (Note that the dashed line in fig 2b is forbidden so this stopping criterion is indeed 
unique). 

2.3.2 The size of 

We want to compute the distribution (|2^). We then split up the phase space integral according 
to the partition into subsets Qq 




(32) 
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where the sum runs over all coprime lattice vectors. We will now smear this distribution 

P-(0 = y"/ Sail ~ Uxs))dxs . (33) 

The smearing functions 6a-{x) may be choosen as gaussians (this is not essential, the important 
thing is that the flanks decay at least exponentially). The mean and variance is given by 



/ xScr{x)dx = 0(a) 
Jx^6a{x)dx = 0ia^) ■ 



(34) 



The exact mean and variance is left undeterminate and may be different from one place to another. 
This will enable us to write equalities such as 6a{x + e) = Sa-{x) as long as e does not exceed a in 
magnitude. 

We will choose a ^ 1. We will eventually consider Pa{l) as a function of ^ = 2RI. The 
smearing width in ^-space is thus ~ R. The goal is to average log^ over p{l{^)). The error 
induced by exchanging p{l) with the smeared version Pa-{1) will be discussed in section 2.4. 

We will also assume some ambiguity concerning the normalization of the smearing functions 

6aix)dx = 1 + 0{R) . (35) 

We will now rewrite eq. (Q) using our conventions on the smearing function. If Xs G Jlq then 
ls{xs) = q + 0{R) and we write 

pAO^Y^f SAi-q + o{R)) = Y.a^6a{i~q) . (36) 

where 

dxg (37) 



and q = |q|. Our goal here is to compute Oq. 

Consider a trajectory hitting disk q. The relation between the phase space variables and a 
and the scattering angle /3q on disk q is (see fig 2). 

sin(0 — 6q — a) + sin a = sin /3q , (38) 

R 

where 9q is the polar angle of the lattice vector q. Provided that there are no interposed disks 
Oq is given by — 1 < sin/3q < 1 and — 1 < sin a < 1. Expanding sin((/) — 9q — a) we get 

^{(l> - 9q - a) + 0{R'^/q'^)+sma = sin f3g , (39) 
R 

It is convenient to change integration variables in dxs = d(l)d(sma) from ((j),sm{a)) to 
(sin(/?g)), sin(Q!)). We thus get 

dXs = — d(sina)d(sin/?q)(l + 0(R^/q^)) . (40) 
47rq 



First consider the case q < 1/2R: 
As we showed in ref. |l0| disk q is then fully accessible and 

a, = ^(l + 0(i?2)) , q<^ . (41) 
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Next consider disks lying beyond the horizon q > 1/2R: 
We make use of the results of section 2.3.1 and write q„ = q' + nqc- We say that disk q„ lies 
in the qc corridor. Note that to each qc the corresponds two q' one with smaller and one with 
larger polar angle (preceeding and succeeding qc in the appropriate Farey sequence). Below we 
assume that q' is the one with the larger polar angle, the other case is completely analogous. The 
geometry in the qc corridor is illustrated in fig 3. If Qq^ were not obscured it would be given by 

%((?!)-a) + 0(i?4)+sina = sin/3q„ , (42) 
it 

with qn = |q' + nqd- Without loss of generality we can set 6^^ = 0. The reason why 0{R'^/qf^) 
in (|39| ) becomes 0{R'^) here is that g„ > 1/27? by assumption. From now on we set /3q„ = /3 as 
there is no risk of confusion. 

The disk q„ under observation is shadowed only by two disks, namely q„_i and qc. The 
question is now where the phase space of these cut through the (a, /3)-plane corresponding to q„. 
First we consider qc. The relevant border is described by an equation obtained from (^9|). 

^(0 + 7„-a) + O(i?4)+sina = l , (43) 
H 



We now combine (42) and ( [4 3D to eliminate (j) and get 

qn{l + 0{R*j) ^ qcSin 13 + (qn - qc) sma+ ^qcqnln ■ (44) 

H 

With some elementary geometry one may show that qcq-nln = 1 + 0{R^) (to get the 0{R^) 
correction one has to use that gc > 1 and qn < l/2i?), so 

(j„(l + 0(i?'^))=(Zc sin/3 +(g„-gc)sina + 1(1 + 0(i?2)) . (45) 

H 

Next we do the same thing with qn-i. The equation for the borderline is 



(0-(7„_i-7„)-a)+sina + O(i?4) = -l . (46) 



R 

Combining (^) and (|2|) we get 

- 9n(l + 0{R^)) = (7„_i(sin/3 - sina) + g„sina - l^n-i'JnlTn-i - In) ■ (47) 

ii 

Another geometry exercise shows that qn-i —qn^Q + 0{R^) and g„_i(7„(7„_i — 7„) = 1 + 0{R^) 
so we arrive at 

-(?„(l + 0(i?*) = (q„ -gc)sin/3 + <7cSina- + ■ (48) 

To leading order in R (^ ) and (^^ describe straight lines cutting through the (sin /? sinQ;)-plane. 
Provided that qc < 1/2R there is a non vanishing area in upper-left corner in the (sin/3, sina)- 
plane for any qn (lower left corner if q' is chosen as the Farey predecessor). So we arrive at the 
following results 

f 1 9 < 2l? 

a, = ^(l + 0(i?2)). Ju<<l<JB+1c (49) 

I (q-q,)(q-2q^) 2R ^ 1c ^ 1 
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2.3.3 The density of coprimes 

The derivation below is very similar to the standard derivation of the asymptotic behaviour of 
sums over Eulers totient function J2n=i for details see e.g. ^. 

We are interested in the number of coprime lattice points inside a certain radius r (in this 
section we do not restrict the consideration to the first octant) 

N{r)= 1 ■ (50) 

m, n 

gcd^m, n) — 1 

We are going to relate this function the total number of lattice points inside the circle of radius r 

i?(r)= E 1 ■ (51) 

m, n 

This is possible since any vector (m, n) may uniquely be written as (m, n) — d ■ {m, fi) where 
d — gcd{m, n) and {rh, n) are coprime. we can thus write 

r 

i?(r)-E^('^/^) ' (52) 

d=l 

Using Moebius inversion theorem we get 

r 

N{r) = E f^{m)R{r/m) , (53) 

m— 1 

where ^i{m) is Moebius function. The hmction R{r) is 

R{r) = Trr^ + E{r) , (54) 

where the error term, according to a a theorem of Sierpinski, is E{r) = 0{r^^^) which means that 
there exist a constant C such that |i?(7')| < C r^^^ for all r > 1. Combining the two previous 
expression we get 



E m(w)— - E li'{Tn)—\^Y^[i{rn)E{rlm) , (55) 

m— 1 m— r+1 / m—1 

where the first infinite sum is related to Riemann's zeta function ^^^i /tti^ — 1/C(2) = 
G/tt^. Now to the second term 

OO - OO ^ 

E — l< E — = Oil/r) , (56) 

m— r+1 m— r+1 

where we used < 1. The third sum may be given a similar bound 

OO oc 1 

E Km)E{r/m)\ < E \E{r/m)\ < C r"'^ E ^ ^ ^ir) . (57) 



m—l 
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So we arrive at the result 

N{r) = -r^ + 0(r) . (58) 

TT 

We will use the derivative of this function, which is the mean density of coprimes with respect to 
the radius 

4(0 = ^-^ + 0(1) . (59) 
ar TT 

One should of course be careful when integrating asymptotic series like this, but the procedure 
above is legitimate for any smearing of the exact N(r). 

2.3.4 Assembling the pieces 

In this section we are going sum up the different contributions computed in section 2.3.2 by 
means of the means density of coprimes computed in section 2.3.3. 

First we consider the case / < 1/2R. 
Then ttq is a function of q only: Oq — a{q) — R/{nq) ■ (1 + 0{R^)), cf. eq. (49) 

p.ii) = E - <?) = 5] - • (1 + o{R^))5,{i - q) 

i — ' ^ — ' Txq 
q q 

= V—j,(;-g)-4(o4 = ^ + ^ -0(1/0 ■ (60) 

q 

We appealed to eq ( p5| ) to get rid of the 0{R^) term in going from the first to the second line 
above. We reformulate the last expression slightly 

pM0)^^ + O{R')-O{1/0 2i?<e<l , (61) 

where ^ ^ 2RI. 

Next we consider the transition region 1/2R < I < 1/ R. 
We then use the result of section 2.3.2 to write q = q' + nqc. According to eq. ( p9[ ) the amplitudes 
Oq depend on the size of corridor qc = jqd and the length of q: flq — a{qc,q)- The length q is 
q = Aqc + nqc + 0{R), where A is defined in fig 3, it is a function of qc. We get 

oo 

p,(/) = 2EEa(g„g)=J,(/-g) (62) 

QC 

= 2EEa(gc,g)(5^(/-gcri-a + 0(i?)) , 

where we inserted the factor 2 to account for both parents q' and q" of q^. The 0{R) in the 
argument gets swallowed by the width of the delta function. Now we will turn the sum over qc 
into an integral over the density of coprimc lattice vectors 

/•1/2R tx> pi 

Pa{l) = 2 dqcd,{q,) a{qc,l)y2 dAfg^{A)S^{l ~ q,n ~ q,A) , (63) 
"'0 . oJO 
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where /^^(A) is the distribution of the parameter A. Since this gets smeared by the width 
a/qc ^ I/qc the results of the Appendixcan be apphed and fq^{^) can be replaced by a uniform 
distribution fq^{A) ^ 1 

I-1/2R i-oo 

Pail) ^2 / dqc dc{qc) a{qc,l) / driS^il - qcV) 

Jo J2 
I.1/2R ^ 

= 2 / dq, d,{q,) a{q,, l)-Ba{l - 2q,) (64) 
Jo Qc 

l.min(l/2R,l/2) 

= 2 dqc dc{qc) a{qc,l)— , 

Jo Qc 

where Oa{x) is a smeared step function. The restriction the the range I < l/R implies that 
■min{l/2R,l/2) = 1/2. Next we insert the expression for a{qc,q) from eqs. (|49|). We then dispose 
of the 0{R^) as before and for the moment suffice with leading term of dc{qc) 

/,^ „ f^'^ , ^'2qc 1 R,, (l/2R-lf , 
Pail) = 2 dqc—^ -(1 - f )+ (65) 

Jl-1/2R Qc TTi qc[l - qc) 



i 



-'/'^ 12qc I R Jl/m-qc)^ 

dqc JY, rr-, TT^ + EiO • 

TT qcTTl (I - qc)[l - 2qc) 



We change integration variable to ry = 2Rqc and use as before ^ = 2RI and we have to solve the 
integral 

24R ( f^/' drj (1-0% ^ f^-' dv jl-rj)^ \ 
(2$ + e(4 - 30 log(0 + 4(e - 1)2 log(^ - 1) - (2 - 0' log(2 - 0) + ^(0 ■ (66) 



Now remains only the case I > l/R: 
The calculation is completely analogous to the previous case and we get 

^l/2R ^ 
P<y{l{0) = 2 1 dqc dc{qc)— a{qc, q = I) 

Jo qc 



12qc 1 R Jl/2R-qc)^ 

dqc JY, rr-i ^ + 67) 

TT qcTTl (I ~ qc){l ~2qc) 

2AR dT] (1 - rj)^ 



6R 



^ (e-^)(e-2ry)+^^^^ 
(2e + e(4 - 30 log(0 + 4(e - 1)2 log(C - 1) - (2 - 0' log(^ - 2)) + E{0 



The 1 < ^ < 2 and ^ > 2 results can be summarized by, 

12R , ^(^) ^ < 1 

(68) 



pMO) -\ + _ log(^) + „ 1)2 iog(^ - 1) - (2 - 0^ log |2 - ^|) + £;(0 ^ > 1 
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where ^ = 2RI. It is instructive to perform a Laurent series expansion of Po-(Z) for ^ > 2 

2*= - 1 1 
fc(fc + l)(fc + 2)^fc+2 



48i? 2*= - 1 1 



fc=i 

We now have an exceUent opportunity to check our calculations by computing the normal- 
ization factor J piy(l)dl. After a minor orgy of integration it is comforting to discover that this 
factor is exactly unity. The expectation value / lpa-{l)dl is found to be 1/2R, which is correct up 
to 0{R). 

Now remains only the error term E{£^). The remaining error is now exclusively due the 0(1) 
term in dc{r) — 12r/7r + 0(1). The case 2R < ^ < 1 is already clear. The case 1 < ^ < 2 is also 
simple since there Oq < R? j^-n so -E(C) — 0{R^)0{£,'^). The remaining case ^ > 2 needs a little 
more care since we would like to keep the 1/^'^ behaviour. 

. drj (1-77)2 ^2Ci?2/-ii_^ 



m)\ < CR' / ' " < — ^ / '-drj = 0{R' log R)0{l/e) ■ (70) 

J2R V ^{^-V){^-2V) J2R V 

We summarize 

r 0{R'^)0{1/0 2i? < ^ < 1 

m)^{ 0(i?')0(l) l<^<2 . (71) 

[ 0{RHogR)0{l/e) 2<<e 

2.4 The small R limit of the Lyapunov exponent 

Now we will finally compute the integral / log£,-p{^)d(^/2R without boring the reader with details. 
When integrate over the region 2 < ^ we use the expansion (^9|) and express the integral in terms 
of the sums 

EZiw- C(2) = -V6 

EZi w = C(3) = 1.2020569031 . . . 

Er=i2^= log 2 . (72) 

T.Z,lh= -Vl2-ilog^2 

Er=i2^= -log2f^ + ilog3 2 + |C(3) 

After this we might as well express the rest of the integral in terms of these sums. The resulting 
c(i?) (as defined in eq (|o|)) is 



c{R)= J\ogipiOd^/2R^ Jlog^p,iOd^/2R + E,{R)^^C{^)-i\og2 + EiiR)+E2{R) , 

(73) 

where Ei{R) is the error induced by smearing the delta functions, cf section 2.3.2, this may be 
shown to be Ei{R) = 0{R\ogR). The second error term is E2{R) = J^\og^E{0d^/2R = 
0(i? log^ R) where the dominating contribution comes from the vicinity of ^ ~ 2R. 
We can now state our main result, cf. eq. (^l|), which is 

A = -21ogi? + C + 0(i?log^i?) , (74) 

with 

C 1 -41og 2 + -^C(3) « -0.12837205 . (75) 
27r^ 

We stress that the derivation in the appendix is still not rigorous and the error estimate 
0(i? log^ R) is only conjectural. 
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3 Discussion 

We stress that the approach the considerations in section 2 appUes to a large class of system but 
the function p(l) need to worked out for each particular geometry. The term —2 logi? appears to 
be universal for two-dimensional Sinai- billiards § |, |. One may use rather crude methods ||] 
to estimate p{l) and get quite decent estimates of the constant C in eq. ([74|). 

The distribution of recurrence times contains a lot of dynamical information, concerning e.g. 
correlation decay and diffusion behaviour |ll|, ^ . 

I would like to thank Predrag Cvitanovic for comments on the manuscript. This work was 
supported by the Swedish Natural Science Research Council (NFR) under contract no. F-FU 
06420-303. 
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Appendix 

In this appendix we will present evidence that the parameter A, as defined in fig 3, is uniformly 
distributed between zero and unity. 

Let q'(q) be the mother (cf. section 2.3.1) of q and q'(q) its father. Further, let A'(q) be the 
A- variable (see fig. 3) associated with the mother and similarly with the father. Let now Fr{A) 
be the (probability) distribution of all A's associated with all coprime lattice points q such that 
g = |q| < r 

N{r) ■ Fr{A) = J2e{r- q)^ {S^{A - A'(q)) + ^„(A - A"(q)) (76) 
q 

where w is a smearing width to be defined later. As in section 2.3.3 we also assume smearing of 
the step function we do not specify it. Let us further define 

d{r)-MA) = -^N{r)-Fr{A) (77) 

where Fr{A) and /r(A) will share their leading asymptotic behaviour. 

Next we use the result of section 2.3.1 to write each q as either q = q'(qc) + nqc or q = 
Ct"(ctc) + nq^c- Let us consider the entire sequence q„ = q'(qc) + «qc- Then q'(qn) = qn-i and 
q"(qn) = qc and 

A'(q„) = + oiwj = ""'^iz:' + ^( V'^") 

|qn| A'(qc) + n 

^"^^"^ = 1^ + ^^'/'^"^ = aWT^ + ""^'^''-^ ^''^ 

Qn = (A'(qe) + n)qc + 0{l/{nql) 

And we can write 

oo 

N{r) ■ Fr{A) = ^ ^ 0(r - (A'(q,) + n)ge + 0(1/ {nql)})- (79) 

^-^^ - + 0{l/{nq,r)) + 6^{A - \ + 0{l/{nq,r)) 

A'(qc)+n A'(qc)+n 

where wc have used the fact that the sequence q = q"(qc) + ?^qc give rise to an equivalent term. 
Next we insert the identities 



J 6{r' - qc)dr' (80) 

/^„(A'-A'(q,))dA' (81) 
Jo 



1 

1 



After some work we arrive at 

N{r) ■ Fr{A) = I dr' I dy9{r - yr')d{r')fr'{y mod 1) • (82) 



((5„(A - ^ + 0(l/r'2)) + (5„(A -- + 0{l/r''))\ 

{ y y ) 
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We now try to solve this integral equation asymptotically. We may restrict ourselves to the case 
A < 1/2 since Pr(A) is symmetric around A = 1/2. We then only need to consider the second 
delta- function in (|8j). We consider now only the contribution from lattice points ri < q < r2 
where ri ~ r2 ^ r is of the same order of magnitude. To get rid of the error term in the argument 
of the delta function we need to take cu = 0{l/r'^) 

N{r2)-Fr,iA)-N{n)-FrA^) (83) 

dr' / dy {e{r2 ~ yr') - e{n - yr')) d(r')/r'(?y mod 1)S^{A - -) 

We insert the leading asymptotic solutions N ~ 6r^/7r and d ^ Ylrj-K and seek the leading 
distribution Fr(A) ~ /r (A) ^ /o(A). The integral equation then gives 



6r| 6r^ 



p(A) (84) 



g^2 g^2\ /•1/2 

^ ^' ' (ix(5^(A - a;)/9(l/a; mod 1) 



which is a smeared version of the functional equation 

p(A) = p(l/A mod 1) (85) 
which combined with the symmetry constarint 

p(l - A) = p(A) (86) 

has the obvious solution p(A) = 1. 

What we have presented is of course not a complete proof but we have presented evidence 
that the distribution of the variable A obtained from all lattice points insider radius r tends 
uniformly to a uniform distribution when r ^ 00 provided that the smearing width uj = 0(l/r^). 
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Figure captions 

Figure 1: a) The Farey tree. b)Part of the Farey three corresponding to the sequence q„ = q'+nqc- 
Figure 2: Definition of phase space variables use in section 2.3.2. 



Figure 3: The geometry of the Qc corridor. 



